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Abstract 

Oscillons are long-lived, localized, oscillatory scalar field configurations. In this 
work we derive a condition for the existence of small-amplitude oscillons (and provide 
solutions) in scalar field theories with non-canonical kinetic terms. While oscillons have 
been studied extensively in the canonical case, this is the first example of oscillons in 
scalar field theories with non-canonical kinetic terms. In particular, we demonstrate 
the existence of oscillons supported solely by the non-canonical kinetic terms, without 
any need for nonlinear terms in the potential. In the small-amplitude limit, we provide 
an explicit condition for their stability in d+1 dimensions against long- wavelength 
perturbations. We show that for d > 3, there exists a long-wavelength instability 
which can lead to radial collapse of small-amplitude oscillons. 
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1 Introduction 



Scalar fields with non-canonical kinetic terms are used ubiquitously in cosmology. They are 
especially prevalent in modeling of the inflaton (e.g. (TJ [2]), dark energy and modifications 
of gravity (e.g. OH]). A significant amount of work has been done on their homogeneous 
evolution in an expanding universe and the evolution of linearized fluctuations about this 
homogeneous background. However, less attention has been paid to their spatially varying, 
nonlinear dynamics (e.g. El [7]). In particular, localized, time- dependent, soliton- 
like solutions in theories with non-canonical kinetic terms have been rarely discussed. If 
such configurations exist, they would be novel objects from a mathematical standpoint. If 
they form a significant component of the energy fraction of the universe, they might have 
cosmological consequences. 

In this paper we show that in a general class of scalar field theories with non-canonical ki- 
netic terms and/or nonlinear potentials, there exist extremely long-lived, spatially localized, 
oscillatory field configurations (oscillonsjB in [HI [10] . While oscillons in scalar field theories 
have been studied extensively in the literature (for e.g. [12l HH HH HS1 [16l HZl HH UHl Efll EH 
E^[231[2S[23[2I3[^ every instance so far uses a nonlinear term in 

the potential and a canonical kinetic term. To the best of our knowledge, this is the first time 
their existence is being demonstrated in the presence of non-canonical kinetic terms. When 
the non-canonical kinetic terms are significant, we refer to the oscillons as "k-oscillons"[j]. 
We derive the condition for their existence, including effects from the non-canonical kinetic 
term as well as the nonlinearity in the potential. Our results are general enough to include 
oscillons supported by nonlinear potentials, oscillons supported purely by the non-canonical 
kinetic terms as well as oscillons supported by a combination of both. The inclusion of 
non-canonical kinetic terms significantly expands the space of theories where oscillons can 
exist. 

Our analysis is done in a small-amplitude approximation, but is otherwise quite general. 
We consider scalar field Lagrangians of the form 



r}^ is a Minkowski space metric with the 'mostly +' signature. The only restriction of T 

* Oscillons are similar to Q-balls [IT] in that their existence has nothing to do with topology, however 
unlike Q-balls, oscillons do not have an exactly conserved charge. 

' Scalar field dark energy with significant non-canonical kinetic terms is often referred to as k-essense (as 
opposed to quintessence). 



T(X,<p)-V(<p), 




in d + 1 space-time dimensions where 



(2) 
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and X is that they can be written as 



T(X,<p) = X + ^ 2 X 2 + ^X 2 + ... 

\r( \ 2 , ^3 3 , ^4 4 , ^5 5 , ^ 

Vw = 2<P +■■■ 

where all the field variables, space-time co-ordinates and coefficients have been made dimen- 
sionless using appropriate scalings. In the next section we motivate this Lagrangian and 
discuss scalings of the parameters and fields in terms of a mass and a cutoff scale. An- 
ticipating a small-amplitude expansion, we have organized the series in terms of powers of 
the field and kept terms up to fourth order in the fields. Note that terms of the form tpX 
and f> 2 X which should be included in the above expression can always be absorbed using a 
field re-definition. Furthermore, note that the X 2 term cannot be eliminated using a field 
re- definition. Our choice of T and V also ensures that we recover a free, canonical scalar 
field theory when X, tp — > 0. This form is general enough to cover a large class of scalar 
field theories of interest in cosmology including axions [35], DBI inflation [Tj, monodromy 
inflation [Ml EI], k-essence [I] and scalar-tensor theories [3]. 

For the Lagrangian discussed above, and in the small-amplitude, spherically symmetric 
case, we provide 

• a condition for the existence of oscillons in terms of a relationship between the first few 
coefficients in the series for T and V, 

• explicit, controlled, analytic solutions in 1 + 1 dimensions and approximate solutions in 
3 + 1 and higher dimensions, 

• a condition for their stability against long wavelength perturbations. The condition 
shows the presence of a long wavelength instability for d > 3 (in the small-amplitude 
approximation) . 

We also calculate the energy loss from these general oscillons due to an expanding back- 
ground. 

For canonical kinetic terms, the authors of [271 [2H Ell [281 I2H] show that oscillons are 
produced copiously at the end of inflation (as well as in phase transitions) and can dominate 
the energy density of the universe at that time (e.g. [2H])- A similar situation is possible 
with a somewhat contrived model of dark energy as well [3H]. Models with non-canonical 
kinetic terms and/or nonlinear potentials are well suited for amplifying fluctuations around 
a homogeneous oscillatory background field. This provides a natural mechanism to amplify 
quantum fluctuations existing at the end of inflation, possibly leading to the formation of 
large-amplitude oscillons. We will pursue this possibility in inflationary models with non- 
canonical kinetic terms in future work. Indeed our original motivation for studying oscillons 
in non-canonical theories arose while trying to explore self-resonance and preheating in the 
DBI scenario (see for e.g. jHHE]) 
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The rest of this paper is organized as follows. In the section 2 we motivate the Lagrangian 
discussed in the introduction, in section 3 we derive an oscillon solution in d + 1 dimensions, 
in section 4 we discuss the stability of the oscillon solutions and in section 5 we summarize 
the main results and discuss directions for future work. 



2 The effective Lagrangian 

In this section we discuss the motivation for the Lagrangian in the introduction. For the 
reader interested in oscillons from a mathematical standpoint, this section can be omitted 
without affecting the rest of the paper. From an effective field theory perspective, a general 
Lagrangian with non-canonical kinetic terms has the form 

X 2 

= a^X* + a 2 (0)-^ + ...- C/(0), (4) 

where = — (l/2)<9£0<9^0. We are assuming that only first derivatives of the field appear 
in the Lagrangian. Let U(<f)) = (l/2)m 2 2 + . . . and A>mis the cutoff scale. Furthermore, 
we assume that a Taylor expansion exists for all a n (<j)) and a\(0) = 1. A field re-definition 
enacted via d(j)/d<j) = a/oi(0) yields 

X- 

where X^ = — (l/2)dp,4>d p, 4 > . We Taylor expand C around (0, X$) = (0,0) and assume that 
as (Xx, 0) — > (0, 0) we recover a canonical free-field Lagrangian. This yields: 

C = Xf + £ 2X 3_ + . . . _ l _ m ^ _ i m ^%0 3 - -m^X^ + . . . (6) 

Let us redefine the fields and space-time variables as follows 

x^ = mx' 1 , 

l-d - 

ip = m 2 0, 



X = m~^X h (7) 



(d+i) 



& = (a J 6. 

With these redefinitions we get the Lagrangian we will use for the rest of the paper: 
C = T(X,<p)-V{<p) 

= [x + 6X 2 + ...]- f^ 2 + ^V + ^V + -- 



(8) 
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where anticipating a small-amplitude expansion, we have only kept terms up to 4th order 
in the field </?. Note that the next order terms in the first and second brackets above would 
have the form £,3<pX 2 and (l/5)X^,(f 5 respectively. 

Another way of arriving at the Lagrangian above is as follows (see for e.g. [12] )• Consider 
a Lagrangian with two fields in 3 + 1 dimensions (the generalization to d + 1 dimensions is 
straightforward) : 

£ = ~d p ^4> - \d-^d^ - U 2 ^ 2 - l -m 2 d? - l -m\,d7 - - \J^\d-^ +... 

(9) 

with A ^> m. The heavy field will sit at the minimum of its effective potential, with the 
corresponding field value given by 

60, 



^uj^y (1Q) 
v V 2 A 3 

Substituting into the original Lagrangian and setting the kinetic term of the heavy field to 
zero, we have 

£ = -\d,m + 6 ( ^? r - \™ 2 <P 2 - - + ... ^ 

With appropriate scalings of space-time and fields with mass m and A defined in equation 
<^j, we recover the effective Lagrangian discussed above and in the introduction (see equation 
©)• 

Note that if £ 2 and A n are order one, we expect £ 2 "C A„ and £ 2 *C 1 if A ^> m. However, 



this need not always be the case. For example, in the DBI case £ = / 1 (<p) 1 — a/1 — 2f(<p)X 

^(^P] yields £ 2 = 1/2. For the rest of the paper we do not make any particular assumptions 
about the sizes of £ 2 and A n , apart from assuming that they are not much larger than unity. 



3 Oscillons in d + 1 dimensions 

Let us begin with the equations of motion associated with the Lagrangian presented in the 
introduction and discussed in the previous section: 

n , a i T a w , dxd v T 8,„V-8 V T 

Op + g^fd^X + sfrO*** = g x T ■ (12) 

We are interested in small-amplitude, radially symmetric, spatially localized, oscillatory (in 
time) solutions. First, we rescale the time and space variables by a small parameter e as 
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follows: 

T = (13) 
p = er. 

and expand the solution as 

<p{t, r) = C 0i(r, p) + e 2 2 (r, p) + e 3 3 (r, p) + . . . (14) 

Plugging the above scalings and form of the solution into the equation of motion ( 1T2|) . and 
collecting terms to lowest order in e we get 



d TT <j)l + 01 = 0, 

=> <j>i( T >P) = /(p)cosr. 
where we assumed that d T 0i(O,p) = 0. At second order in e we get 

d T r4>2 + 02 = ~A 3 / 2 (p) [1 + COS 2t] , 

2 (r,p) = ^A 3 / 2 (p)[-3 + 2cosr + cos2r]. 



where we assumed that 2 (O,p) = d T 4> 2 (0,p) = 0. At the next order in e, we get 

<9 TT 03 + 03 



(15) 



(16) 



dlf + ^d p f-f + l^ 2 -x i + ^xl)f 

+ [. . .]cos2r 
+ [■■■] cos 3t. 



COST 



(17) 



If the coefficient of the cost term is non-zero, 03 would grow linearly with time, inconsistent 
with the time-periodic solution we are looking for. @ To avoid linear resonance, we need to 
set the coefficient of cost to zero. This in turn yields the equation for the profile f(p): 



d 2 p f( P ) + ^-^-d p f{p) - f(p) + ^A/ 3 (p) = 0. (1* 



where for future convenience we have defined 

A = 6 - A 4 + yA 2 , (19) 

whose sign will turn out to determine whether oscillons exist or not. 

This profile equation ( IT8|) is valid in any dimension d, but can be solved exactly for d = 1. 
Let us consider the d = 1 case first. 



*Note that the [. . .] cos2r and [. . .] cos3r terms will yield a periodic solution for 03 . This solution can 
then be used in calculating terms at the next order, just as the periodic solution for (f>2 played a role in the 
03 equation. This pattern extends to all orders. 
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3.1 1+1 dimensional oscillons 



For d = 1 the profile equation becomes 

3 
4 

One can think of the above equation as describing the motion of a particle in a potential 



d 2 p f(p)-f(p) + -Af(p)=0. (20) 



U(f) = —— + —A/ 4 . (21) 
2 16 

The energy associated with this motion is conserved, and given by 

S=U(f )= { -^+U(f), (22) 

where we have used d p f(0) = and /(0) = /q. Now since the solutions are localized, the 
energy £ — U(fo) — 0. This immediately, yields 



For a localized solution to exist, we need 

Tn 

(24) 



A = 6 - A 4 + ^A 2 3 > 0. 



This is one of our main results. Before moving on to solving the profile equation, we pause 
to discuss A in a bit more detail. If the non-canonical terms are absent (£2 = 0), we need the 
usual "opening up of the potential" condition: — A4 + (10/9) A3 > to get oscillons. More 
importantly, note that for a quadratic potential (i.e. A n = 0), the non-canonical terms are 
sufficient to yield oscillons. For example if A n = 0, then £ 2 > is sufficient. It is also worth 
noting that in models with non-canonical kinetic terms, the sound speed differs from 1. For 
the model under consideration, the sound speed is 

c 2 = (l + 2X d2xT 



d x Tj (25) 
= 1 - 4£ 2 X + . . . 

Thus, for A n = 0, the condition for having k-oscillons is the same as the sound speed being 
less than 1. 

Now, let us get back to solving the equation for the profile f(p). Using /(0) = f derived 
above and integrating (d p f) 2 /2 + U(f) = yields 



/(p) = \/^sech(p). (26) 
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Explicitly in terms of the original variables: 



<p(t,r) 



— — sech (er) cos 

•J Z_A 



(Vl^t) + 0\e 2 



(27) 



where A = £2 — A4 + yA§. This is our second main result: a small- amplitude oscillon in 
1 + 1 dimensions in theories with non-canonical kinetic terms. This solution has the same 
functional form as the canonical small-amplitude oscillon, apart from the appearance of £2 
in A. One can go beyond the leading order as well: 



<p(t, 



3A 



sechp cos r + e 2 -^^-sech 2 p [—3 + 2 cos r + cos 2r] 



(28 



+ 0[e 



where p = er and r = yl — e 2 t. Note that at second order, only the odd term contributes. 
If the Lagrangian had ip — > —ip symmetry, the the correction is higher order. 



3.2 3+1 dimensional oscillons 

For d = 3, the profile equation becomes 

d 2 p f{p) + - p d p f(p) - f(p) + \M\p) = 0. (29) 

We can view the above equation as describing the motion of particle in the presence of 
a potential U(f) as in the 1 + 1 D case, but now we have a 'friction term' of the form 
(2/p)dpf(p). As a result, the energy 



is no longer conserved. It changes with p as 



d p s( P ) = - 2 -{d p ff. (31) 

With the requirement that the solution is localized, we need £ — > as p — > 00. Requiring 
that the solution is smooth at p = implies d p f(0) = 0. This implies that for a localized 
solution we must have £(0) = U(fo) > 0. Numerically, one finds a localized solution for[§ 

fo - M. (33) 



^more precisely 



fo w 3.06699 x J— (32) 
V 3A 
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Thus the solution will have the form: 

<p(t, r) « ef F(er) cos (y/l - e 2 t) + 0[e 2 ] (34) 

The profile F(p) looks 'sech' like with F(0) = 1. An excellent approximation (at the few % 
level in the 'core' region, deteriorating in the tails) is given by 



F(p) = ^sech (/J (2A)1/3 p) 1 < A < few (35) 
An identical analysis can be carried out in d ^ 1,3 and we do not repeat the derivation here. 

4 Stability 

We now turn to the question of stability of the oscillons against perturbations. To this end, 
we linearize the equation of motion (|T2|) around the small-amplitude oscillons solutions <p osc 
as follows: 

50 - G&VoscVoscStp + [- (l - 2£ 2 0L) V 2 + 1 + 2A 3 y? osc + 3\ 4 ip 2 osc - 3£ 2 y?L] Sl P = °- ( 36 ) 

We have used the following to arrive at the above equation: 

0[ip osc ] ~0[<p 
0[\V<p oac \] ~ e 



(37) 



and kept terms up to order e 2 . The ^ 2 < ^ 2 S c^ 2 ^V 9 term has to be kept since we have not (yet) 
restricted ourselves to long wavelength perturbations. Let us remove the term with the linear 
derivative by redefining the 5ip as follows: 



6<p = xexp 



3£ 2 / dsip 







(3J 



With this redefinition, the equation of motion becomes 

X +[-(!- 260L C ) V 2 + 1 + 2A 3 ^ osc - 3(6 - \M SC ] X = 0. (39) 

The solutions (f osc are periodic in time. Hence stability can be determined via a Floquet 
analysis. It is tempting to Fourier transform the above equation and try to determine the 
Floquet exponents (growth-rate) mode by mode. However, since the background varies in 
space, the Fourier modes do not decouple. A stability analysis in Fourier space, while possible 
(see [26]). is numerically intensive, especially in higher dimensions. We do not follow this 
approach here. Instead we will carry out a stability analysis in position space, but restrict 
ourselves to perturbations which vary on length scales comparable to the size of the oscillon. 
Even without considering short wavelengths, we will show that there exists an important 
instability in dimensions > 3. 



S 



4.1 Long- wavelength stability analysis 

When considering perturbations with wavelengths comparable to the size of the oscillon, we 
can drop the £,2<Posc V 2 5y? ~ f4 term in equation (1391) since it is higher order in e. This yields 



At this point, the equation of motion for linearized fluctuations is identical to that in the 
canonical case, apart from the coefficient £2- From now onwards, our long wavelength stabil- 
ity calculation closely follows the one by Amin and Shirokoff [25], where we related oscillon 
stability to the stability criterion derived by Vakhitov and Kolokolov jH] in the context of 
light focusing in a nonlinear medium. Here, apart from showing this relationship in the con- 
text of non-canonical oscillons, we also provide a pedagogical proof of the stability criterion 
itself in an Appendix (not provided in [25] )• 

Recall that the oscillon solutions have the form 

<Poac(t,r) = ef(p) cos(r) 

+ e 2 ^/ 2 (p) [-3 + 2cosr + cos2r] + ... 
o 

where p = er, r = yl — e 2 t and / is the radial oscillon profile satisfying equation ( ITS]) . 

We are interested in determining the stability of the above oscillons to perturbations 
with wavelengths comparable to the size of the oscillons. With this in mind, let us define a 
scaled spatial co-ordinate p = er where e = ej^fa with a being an order 1 parameter. We 
expect the most unstable perturbations to oscillate at the oscillon frequency: y/1 — e 2 with a 
'slowly- varying', time-dependent envelope driven by the oscillating background. To capture 
the time dependence of the envelope we define a slow time T = e 2 t which is to be treated 
independently from r. Note that although 0[e] ~ 0[e], they are independent, e plays the 
role of determining the oscillon solution whereas e is used for analyzing the stability about 
this solution. This distinction is made explicit via the introduction of the a parameter. With 
these definitions we have 



X + [-V 2 + 1 + 2\ 3 <p osc - 3(£ 2 - A 4 VL] X = 0. 



(40) 



dt 2 



(1 - e*)dl + ~e 



e 2 VT^e 2 2d T d T + e 4 <9 2 



V 2 




(42) 



Furthermore, let us expand the perturbation \ i n powers of e as follows: 



X = Xo + exi + e 2 X2 + • • • 



1 1 2 

= Xo + —^tXi + -e X2 + ■ ■ ■ 



(43) 
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We now substitute the space-time scalings, the e expansion of the perturbation and the 
oscillon solution into the EOM of the perturbation ( I40p . Collecting terms order by order in 
e we have (up to order e 2 ) 

[d 2 + 1] xo = 0, 

[d 2 T + l] Xi + [2A 3 $cosr] X o = 0, 
[d 2 T + l] x 2 + [2A 3 $cosr] xi 



28 T d T - V 2 - ad 2 



£2 _ A 4 



2 
3' 



h ^A 2 $ 2 COS T - 



6 - A 4 - --3 



A 2 $ 2 cos2r 



where we have defined 



$(a,p) = y/af(y/ap) = y/a f(p). 



Xo = 0, 
(44) 

(45) 



Note that treating 0[V 2 Xo] ~ ^[Xo] we are restricting ourselves to perturbations that vary 
spatially on the scale of e. The first equation of (jUJ) yields 

(46) 



Xo = u(T, p) cos r + v(T, p) sin r. 
Substituting xo into the second equation and solving for xi we get 



1 



Xi 



A 3 $ [u(-3 + cos 2r) + u sin 2r] 



(47) 



where we have ignored the homogeneous solution of Xi- Finally, substituting xo an d Xi into 
the third equation we get 



[dl + 1] X2 



2d T v 



V 2 - a + ^A$ 2 ) 



COST 



-2dxu ■ 
+ [...] cos 3r + [. . .] sin 3r 



V 2 - a + -A$ 2 ) 



sin r 



(48) 



where A = £2 — A4 + (10/9) A§ (the combination which appears in the oscillon solution). 
Avoiding secular growth requires 

dru = Hiv, , , 

49 

a T v = —H 2 u, 



where 



- V 2 -a + -A$ : 



Hi 



H 2 = --[ V 2 -a + -A$' 



(50) 
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We are interested in the eigenvalues of the above linear system. Let v(T, p) = e v e (p) 
and u(y,T) = e nT u e (p) where (u e ,v e ) is an eigenvector of equation (l4T)j) . Substituting 
(u,v) = e nT (u e ,v e ) into our linear system, we get 



Vtu e = HiV e , 
Qv e = —H 2 u e . 



Equivalently 



-Vt 2 u e 



H 2 H lVe . 



(51) 



(52) 



Since Hi and H 2 are real operators, the eigenvalues — Q 2 are real, that is Q is purely real or 
imaginary. There exist exponentially growing modes if and only if min[— Q 2 } < 0. Hence, we 
now try to determine min[— Q 2 ]. 

In the above form, our problem becomes similar to that of light focusing in a nonlinear 
medium as analyzed by Vakhitov and Kolokolov jH]. Following their techniques, we will 
show in the Appendix that 



where 



sign [min[— f2 2 ]] = sign 



N = / <5> 2 (a,p)d d p 



dN 
da 



(53) 



a J f(V^p)d d p 

a (l-d/2) f f( p)d d p 



(54) 



(l-d/2) 



oc a 

The proof of the relationship between the sign of ^ and sign [min (-fi 2 )] that we used in 
the first line above is somewhat involved, which is why we have moved it to an Appendix. 
Here, we discuss the important and interesting consequences of the result. 

For d > 2, we will have dN/da < and thus sign [min[— Q 2 ]] < 0. Whereas for d < 2 we 
have dN/ da > and thus sign [min [-fi 2 ]] > 0. Evidently, our oscillons are stable against 
long wavelength perturbations in d = 1,2 but not so in d > 2. This is confirmed by our 
numerical simulations. For £ 2 = A3 = 0, and in the small- amplitude limit, this reduces to 
the result in 1251. 



Recall that a is merely a scaling of e in the oscillon solution. In terms of e, the stability 
condition is as follows: Oscillons are stable if and only if 



dN 

d,e 



>0, 



(55) 
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where 



N = 



e 2 Jf(er)d d r = e 2 - d Jf(p)d d , 



p oc e 



d-2 



(56) 



where / is the oscillon profile, that is <p, 



osc 



ef(er) cos y/1 — e 2 t. We have numerically 



verified the existence of a long wavelength instability for d > 3. 

We stress that the stability criterion (1551) is applicable in the small-amplitude limit. More 
precisely it is applicable when a single frequency solution is a good approximation to the 
true solution. We now connect the above result to some related work in the literature. 

In the example of [25] with canonical kinetic terms, unlike our discussion here, the coef- 
ficient of the y? 6 term was assumed to be unusually large. This allowed for an (approximate) 
single frequency solution for the entire allowed amplitude range, which in turn allowed for 
the derivation of the same stability criterion derived above. However, unlike the above case, 
in the large y? 6 case, N was a non-monotonic function of e in 3 + 1 dimensions, allowing 
stable solutions to exist at large amplitudes. We also note that a similar stability criterion 
in term of the oscillon energy was also conjectured in [22] based on numerical results in a 
massless dilaton + scalar field oscillon. While no general stability condition exists for the 
general, large-amplitude case, stability and lifetime of large amplitude oscillons is often in- 
vestigated using Gaussian initial profiles with varying widths and amplitudes. For a flavor 
of such investigations see for example [23l [30] . 

4.2 Radiating Tails in an Expanding Universe 

While we have discussed the stability of our oscillons against 'external' long-wavelength 
perturbations, even without external perturbations, oscillons are not exactly stable. Similar 
to the canonical case, our more general oscillons possess a radiating tail which we expect 
to be highly suppressed [151 I2~T] . with a decay rate ~ e _1//e . Nevertheless, in an expanding 
universe, this tail can be significantly enhanced [1HJ 032 [22] • We briefly sketch out the 
energy loss due to this radiating tail in an expanding universe below. 

For simplicity we will only consider the case in 1 + 1 dimensions. In local co-ordinates, 
the metric for a homogeneous and isotropic expanding space can be written as (space and 
time are measured in units of m _1 ) 



where we assume that H = constant and 0[H/m] = 0[e 2 ]. In this case the solution takes 
the following form (following the technique used in [335]): 




(57) 





(58) 



quantum treatment of the radiation will also increase the decay rate, with the decay rate becoming a 
power law in e [26] . 
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and 



<p(x, t) « e 3 / 2 e -w2y^^cos (sfl - eH - ^x 2 H^j ± < x < JL (59) 
which leads to an energy loss (averaged over time) given by 

— *> ' ¥ « 1 « 77' (60) 



Our analysis appropriately generalizes the result of [19] . The A appearing above contains £2 
from the non-canonical kinetic term along with A3 and A4 whereas in [19], A = — A4. Note 
that this analysis in only valid when m^> H . The energy loss while enhanced compared to 
the Minkowski case, can still lead to lifetimes ^> H . 



5 Discussion 



In this paper we have shown that oscillons can exist in a significantly larger class of scalar 
field theories than previously shown. For a rather general class of scalar field Lagrangians 
of the form ([T]) and 03), we have demonstrated the following: 

• For small-amplitude oscillons to exist, A = £ 2 — A4 + (10/9)A| > where £2 is the 
coefficient of the non-canonical part of the kinetic term. 

• The oscillon solutions in d + 1 dimensions have the form cp(t, r) = ef(er) cos y/1 — eH + 
0[e 2 ] where f(er) is the radial profile. In 1 + 1 dimensions, f(er) = A/8/(3A)sech(er). 
We also provided an approximate form for f(er) in 3 + 1 dimensions. These solutions are 
identical to those in theories with canonical kinetic terms, apart from the appearance 
of £2 m A. 

• The solutions are stable against long wavelength perturbations if and only if dN / de > 0, 
where N = e 2 ~ d J f 2 (p)d d p. 

We have also calculated the energy loss from oscillons due to an expanding background. 

There are a number of natural extensions of our results. The stability criterion above 
is related to long-wavelength instabilities, which we believe to be the most dangerous in- 
stabilities. However, as discussed in the stability section, a calculation of the 'Floquet' 
instability rates at shorter wavelengths while numerically intensive, is also possible. A fur- 
ther detailed investigation of the suppressed radiating tail, effects of expansion on lifetimes, 
'Floquet instabilities' as well as a quantum treatment for these non-canonical oscillons would 
be interesting. 

We have concentrated on the small-amplitude regime in this paper. However, large- 
amplitude oscillons (with large energies) in 3 + 1 dimensions are interesting due to their 
possible relevance in cosmology (see for e.g [271 I2H ESI ESI ES])- In addition, in 3 + 1 
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dimensions, single field, small-amplitude oscillons can collapse due to perturbations with 
wavelengths comparable to the size of the oscillons. As we move to larger amplitudes this 
instability can disappear (e.g. [23]) • Thus an analysis of the large-amplitude case is certainly 
worth pursuing. However, moving to large amplitudes also requires a larger number of 
terms in T and V, which in turn requires a case by case analysis of the solutions and 
their stability. It is of course possible to analyze them numerically. Although we have not 
presented the results here, we have analysed large-amplitude oscillons in the DBI Lagrangian: 

C = 1 — a/1 — 2f(ip)X — (f 2 /2, with rather intriguing dynamics appearing at large 

amplitudes 47]. To keep our analysis as general as possible, and in an effort to present 
analytic rather than numerical results, we have restricted ourselves to the small-amplitude 
case in this paper. An analysis of large-amplitude k-oscillons and their implications in a 
cosmological context is in progress [4"Tj . 
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Appendix 



Proof of sign[min[— Q 2 ] = sign [d N/ da] 

Our proof will closely follow the stability analysis of Vakhitov and Kolokolov presented in 
the context of light focusing in a nonlinear medium jH]. We will need two technical results 
regarding the Hermitian operators Hi and H 2 . 

(a) (u e \H± 1 \u e ) is positive definite for f2 ^ 0. U 

(b) H 2 has only one bounded eigenmode with a negative eigenvalue, all other eigenvalues for 
radially symmetric eigenmides are greater than zero, and the lowest angular eigenmode 
has zero eigenvalue. 

We will assume these to be true for the moment and proceed with the proof. After the proof, 
we justify (a), but for (b) we refer the reader to 



We believe that some of the statements in the rest of the proof can be understood more 
readily based on our experience and intuition with single particle quantum mechanics. In 
particular, H\ (and H 2 ) can be thought of as a non-relativistic Hamiltonian of a particle in 
a finite, radially symmetric potential in d spatial dimensions. Hence we use language from 
quantum mechanics where appropriate. The energies and eigenstates of Hi are denoted by 
{E/3, ^p} and those of H 2 by {£ 7 , 

From equation (1521) we get — Q 2 (Q\u e ) = (Hi<&\H 2 u e ) = 0. The first equality uses the 
Hermitian nature of Hi whereas the second follows because Hi$ = is simply the oscillon 
profile equation (TT8l) . Hence for fl 7^ 0, (&\u e ) = 0. If (M e |if ] ~ 1 |u e ) is nonzero, we can rewrite 
the first equation in (152"]) as 



^ 2 = "e/ ^ ($| ne ) =0 . (61) 





\H 2 \ 


Ue) 


{u e \ 


Hi 1 


\u e ) 



Now, since (u e \Hi l \u e ) is positive definite based on (a) stated above, we have 

sign[min[— f2 2 ]] = sign[min [(tt e |iJ 2 |We)]], (62) 

with 

($|u e ) = and (u e \u e ) = 1. (63) 

We introduce Lagrange multipliers S and /3 to minimize {u e \H 2 \u e ) subject to the above 
constraints: 

F[u e ,S,P} = (u e \H 2 \u e )+S((u e \u e )-l) + P($\u e ). (64) 
The extremum of T is obtained if u e satisfies 

H 2 u e = £u e + (65) 



"where (...)= J ■ ■ ■ d d p and we are using the usual bra-ket notation 
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Moreover, the minimum value of (u e \H 2 \u e ) is given by the smallest eigenvalue of H 2 consis- 
tent with (3>|w e ) = 0. Let £ m ; n denote this minimum eigenvalue. Then 

sign[min[-f2 2 ]] = signfmin [(u e \H 2 \u e )]] = sign[£ min ]- (66) 

We will now try to determine the sign of £ m ; n . As mentioned earlier, one can also think 
of H 2 as the Hamiltonian of a non-relativistic particle in a finite potential well V 2 (p) = 
(a/2) [l — (9/4)sech 2 (y / ap)] . Let {^7} be the eigenstates of H 2 with energies £ 7 . Let us 
expand u e and $ in terms of these eigenstates as $ = ^ 7=0 a 7 ,? / ; 7 an d u e = X^o^y^V 
Plugging these into fl6^]) and using (&\u e ) = we get 



7=0 -T 

where a 7 = (^7!$) and we have defined 
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f3j2f^-: = f39(£) = 0, ( 6 7) 



g{£) = Y.Ph- (68) 



7=0 



Now either (3 = or g(£) = 0. If /3 = 0, then from equation ( 1651) we see that min[£] = £ min is 
obtained if u e = i/)q: the ground state of H 2 , which is radially symmetric and has no nodes. 
This contradicts ($|w e ) = 0. Hence /3 7^ and we have g(£) = 0. We need to find the 
smallest root of g(£) = 

We will now make use of the technical properties of H 2 specified at the beginning of 
this appendix to analyze the minimum value of £ that satisfies g(£) = 0. For the lowest 
radially symmetric eigenstate (ground state of H 2 without any nodes), a = (^ol^) 7^ 0. 
Moreover, from (b), £q < 0. For any radially asymmetric (angular) eigenstate a 7 = since 
that eigenstate will be orthogonal to <3>. In particular if ipi is the lowest angular eigenstate, 
then ai = (^ll^fr) = and by our assumption (b), has £\ = 0. If tp 2 is the next radial 
eigenstate, then by (b), £ 2 > 0. Now consider the behavior of g(£) for £ < £q. In this 
domain, g(£) > 0. For £q < £ < £ 2 , g(£) varies monotonically from —00 to +00 and crosses 
for the 'first' time. Hence if g(£) = in this domain, the root £ is the smallest root £ m mEl] 
Moreover, since g(£) varies monotonically from —00 to +00 in this domain, the sign of g(0) 
determines the sign of £ m i n - Explicitly g(0) > -<=>- £ min < 0. Hence from fl62|) 

sign[min[-fi 2 ]] = sign[£ min ] = sign[-#(0)]. (69) 

Finally, let us now relate g(0) to d(§\§)/da as follows: 

da 

d& /™x 
H 2 — + $ = 0, (70) 
da 



da 



* g{£) varies monotonically between every consecutive pair of distinct £ 7 . 
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Multiplying both sides by $ and integrating, we get 

ld<$|$) 



2 da 



E 

7=0 



(71) 



= -9(0). 

Thus using the above result and equation ( 169|) . we finally have 



sign[min[— Q 2 ]] = sign 


'dN' 

da 


where iV = ($ $) = j 





(72) 



Let us now turn to the justification of the property (a) of Hi assumed in the proof. We will 
show that (u e \Hi \u e ) is positive definite. Note that the eigenvalue problem Hi^/ p = Ep^/p 
is the time independent Schrodinger equation for a particle of mass m = 1 in a radial 
potential well V\(p) = (a/2) [l — (3/4)sech 2 (A/ap)] . Using the profile equation ( [18]) we get 
Hi$> = where $ has no nodes. This implies that \l/o = $ is the unique ground state of 
H\ (up to a normalization) with energy E$ = and all other eigenvalues must be greater 
than 0. Moreover one has the orthonormal set of excited states {^p} with (3^0 which 
satisfies {typ\<&) = 0. For any state which belongs to this subspace spanned by typ, the 
operator Hi is positive definite. Hence iff 1 exists on this subspace and is also positive 
definite. This follows from (^fp\H^ |\& 7 ) = E^ l 5p 1 . Now, from equation fl52|) . note that 
— 2 ($|w e ) = {Hi^\H 2 u e ) = 0. Hence for fl ^ 0, u e lies in the space spanned by {^p}. 
Thus (u e \Hi 1 \u e ) is positive definite. We still need to show that H 2 has only one bounded 
eigenmode with a negative eigenvalue, and the lowest angular eigenfunction has eigenvalue. 
This is somewhat involved, and we refer the reader to [H] where this is discussed further. 
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